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Abstract 

Improving the frequency precision by synchronizing a lattice of oscihators is studied in the phase 
reduction hmit. For the most commonly studied case of purely dissipative phase coupling (the 
Kuramoto model) I confirm that the frequency precision of A'^ entrained oscillators perturbed by 
independent noise sources is improved by a factor N as expected from simple averaging arguments. 
In the presence of reactive coupling, such as will typically be the case for non-dissipatively coupled 
oscillators based on high-Q resonators, the synchronized state consists of target like waves radiating 
from a local source which is a region of higher frequency oscillators. In this state all the oscillators 
evolve with the same frequency, however I show that the improvement of the frequency precision is 
independent of N for large N, but instead depends on the disorder and reflects the dependence of 
the frequency of the synchronized state on just those oscillators in the source region of the waves. 

PACS numbers: 05.45.Xt, 87.19.lni, 89.75.Kd 



I. INTRODUCTION 



Oscillators - devices producing a periodic signal at a frequency determined by the char- 
acteristics of the device and not by an external clock - play a crucial role in much of modern 
technology, for example in timekeeping (quartz crystal watches), communication (frequency 
references for mixing down radio frequency signals), and sensors. A key characteristic is 
the intrinsic frequency precision of the device, which can be quantified in terms of the line 
width of the oscillating signal, or in more detail by the spectral density of the signal in the 
frequency domain, or by the Allan deviation in the time domain. This is the fundamental 
issue, broadly common to all oscillators, considered in the present paper. There are other 
important practical characteristics, including the robustness of the frequency to environmen- 
tal perturbations such as vibrations and temperature fluctuations, that are more dependent 
on the details of the device implementation; these are not considered here. 

Unlike a resonator driven by an external oscillating signal, where the line width of the 
spectral response is determined by the dissipation in the resonator, the line width of an 
oscillator is only nonzero in the presence of noise. An oscillator is mathematically described 
by a limit cycle in the phase space of dynamical variables, and the line width of the signal 
corresponding to a limit cycle is zero. Dissipation serves to relax the system to the limit 
cycle, which is itself determined by a balance of energy injection and dissipation, and does 
not broaden the spectral line. The spectral line is broadened only if there is some stochastic 
influence that causes the phase space trajectory to fluctuate away from the limit cycle. Thus 
the frequency imprecision (line width) is due to noise [1]. 

One way to improve the frequency precision of oscillators that has been suggested in 
various scientific disciplines [2], [3], is to sum the signal from a number N similar oscillators: in 
this case, if the noise sources are uncorrelated over the individual oscillators, simple averaging 
suggests that the effective noise intensity will be reduced by a factor Of course, due 

to fabrication imperfection the isolated oscillators will have slightly different frequencies, 
and so simply averaging the summed signal from the individual oscillators will also tend 
to average the signal to zero. In addition, the line width of the reduced intensity signal 
will reflect the frequency dispersion of the devices. If however some coupling is introduced 
between the oscillators, they may become synchronized to a state in which all the oscillators 
are entrained to run at a single common frequency!^. In this case, the averaging argument 
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would suggest a reduction in the effective noise, and so a factor of N enhancement in 
the frequency precision. 

In this paper I investigate the improvement of frequency precision due to synchronization 
in canonical models of the phenomenon. I focus on lattices of oscillators with nearest neigh- 
bor coupling. I confirm that the factor improvement applies exactly in the case of purely 
dissipative coupling between the phases of the oscillators. However, this scaling breaks down 
if there is in addition a reactive (non-dissipative) component of the coupling. Examples of 
reactive coupling are a displacement, rather than velocity, coupling for arrays of mechanical 
oscillators OE] or trapped ions[7]. In this case the improvement factor becomes independent 
of for large A^, but depends on the amount of disorder. This may lead to much poorer 
frequency precision than anticipated from the naive averaging argument. 

The main focus of the paper is on d-dimensional lattices of oscillators with nearest neigh- 
bor coupling. I briefly discuss the extension to longer range coupling and to complex 
networks. These results should be relevant to questions of the precision of synchronized 
oscillations in biological contexts. 

II. MODEL 

The model I consider is nearest neighbor coupled phase oscillators |H I5HTT] 

ei = u, + Y, r(^j -0^) + Ut), t = i,...,N. (1) 

nn j 

Here Ui are the frequencies of the individual oscillators, which are assumed to be independent 
random variables taken from a distribution g{uj) with width ct[12J. The term ^j(t) represents 
the noise acting on the ith oscillator. I will assume white noise, but the results are easily 
generalized to colored noise, and also to noise that depends on the phase 6i. An important 
assumption is that the noise is uncorrelated between different oscillators. Thus 

m)ut')) = cit-m,, (2) 

where for white noise 

c{t-t') = f5it-t'), (3) 

with / the individual oscillator noise strength, taken to be the same for all oscillators, since 
I am imagining a system where the oscillators are designed to be as similar as possible. 
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The coupling between the oscillators is given by the function F, a 27r periodic function of 
the phase differences of nearest neighbor oscillators. A commonly used model[Hl E] is given 
by the coupling function 

r(0) = sin0. (4) 

Any parameter K multiplying sin0 and giving the strength of the coupling may be scaled 
to unity by rescaling time and frequencies. Thus the only parameters defining the behavior 
of the system are the distribution of the frequencies cuj, and in particular the width a of the 
distribution after this rescaling (i.e., the width of the frequency distribution relative to the 
coupling strength). 

The coupling function Eq. Q is antisymmetric r(— 0) = — r(0), and equation ([T]) is 
purely dissipative[l3]. A more general couphng function 

r(0) = sin0 + 7(1 — COS0) (5) 

breaks this symmetry for nonzero 7, and includes non-dissipative, propagating effects [TOt HI]. 
I will use this model to study the effect of reactive coupling. Without loss of generality, I 
take 7 > 0. 

For small disorder, the phase difference between nearest neighbor oscillators will be small. 
A convenient approximation for F good for small phase differences is [13] 

F(0)^7-'(e^'^-l), |0|,|70l«l. (6) 



III. SYNCHRONIZATION 

I first describe the behavior predicted by Eqs. ([l]|5]) in the absence of noise. For sufficiently 
weak disorder (small a) and for a finite number of oscillators, the oscillations described by 
Eqs. (1-3) with/ = become entrained in the sense that all the phases advance at the same 
constant rate 

Oi = n. (7) 

The solution to these equations can be written 
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(8) 



with 6l a fixed point solution in the rotating frame, and the phase of the collective limit 
cycle given by 



when solutions exist. The behavior in the thermodynamic limit N ^ oo for different 
lattice dimensions d, and the critical value of cr for the onset of this entrained state and 
its dependence on system size, lattice dimension, frequency distribution etc., are subtle 
questions that have not been fully answered. But for the practical case of a finite number of 
oscillators, we expect that an entrained (fully frequency locked) state will exist for sufficiently 
small (t[T6]. Such a state is a limit cycle of the system of oscillators with frequency Q. By 
moving to a rotating frame, 6i{t) — t- 6i{t) — Qt, the entrained state becomes a fixed point, 
simplifying the subsequent analysis. 

The nature of the synchronized state depends sensitively on whether the coupling is purely 
dissipative (7 = 0) or also contains a reactive component (7 7^ 0). For purely dissipative 
coupling the interactions cancel when summed over a block of oscillators. This means that 
the frequency of the entrained state is the mean u of the oscillator frequencies. For a one- 
dimensional lattice, the individual phases 6^'^'^ are then given by 6^^^ — 6^'^'^ = — sin~^ Xi 
with Xi = Yl]=ii^j ~ The accumulated randomness Xj performs a random walk as 

a function of the lattice index i, and so the strain 6',-^^ — 6^^^^ also varies with i roughly as 
a random walk. The break down of the synchronized state as the disorder or system size 
increases occurs when the excursion of Xi exceeds unity (remember the coupling strength 
is scaled to one). This occurs for a ~ X^^/^. On the other hand, for 7 7^ 0, the interaction 
terms summed over a block of oscillators do not cancel, and the same arguments cannot be 
made. For 7 > it is found that the entrained state takes the from of quasi-regular waves of 
some average wave length A propagating away from a unique source in the system, located at 
a cluster of higher frequency oscillators [HI [15]. The derivation of this result is described in 
more detail below. In this wave state the frequencies of all the oscillators remain entrained, 
even though the phases vary by more than 27r over the system for X < N, and by many 
factors of 27r for X N. In two dimensional lattices, roughly circular "target" waves are 



e(t) = nt + Q, 



(9) 



with 00 an arbitrary constant. 6i and are given by solving 




(10) 
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found for 7 7^ 0, with the waves propagating away from a source with location again given 
by a core region of higher frequency oscillators [HI [T3] . 



IV. FREQUENCY PRECISION AND NOISE 



The noise terms in Eqs. ([!]) will lead to deviations of the solution from the limit cycle 
Eqs. (8]|9 ) and so to a broadening of the spectral hues of the output signal from the oscillator. 



The full noise spectrum depends on a complete solution of Eqs. ([T]). However, for frequency 
offsets from the no-noise peaks in the power spectrum that are small compared with the 
relaxation rates onto the limit cycle, the effects of the noise can be reduced to a single 
stochastic equation for the limit cycle phase G that gives the collective behavior of the 
entrained oscillators. This result has been derived for a general limit cycle by a number 
of authors using a variety of formalisms [H IT8ti20] . The key idea is that a change in 
corresponds to a time translation, and so gives an equally good limit cycle solution: thus a 
perturbation to O does not decay, and this represents a zero-eigenvalue mode of the linear 
stability analysis[21j. The remaining eigenvalues of the stability analysis will be negative, 
corresponding to exponential decay onto the limit cycle. For time scales longer than these 
relaxation times, it is only the projection of the noise along the zero eigenvalue eigenvector, 
that is important: the other fluctuation components will have decayed away. 

For the white noise sources considered here, this stochastic equation for the phase is 
simply pUj 

e(t) = n + s(t), (11) 

with 

{E{t)E{t'))=F6{t-t'), (12) 
with F the noise strength resulting from the projection and Q the limit cycle frequency 



which is Q with an 0{F'^) correction. The solution to Eq. (11) is a drift of the mean phase 
at the rate 

(e(t)) = Ot, (13) 

together with phase diffusion 

((e(t) - Qtf) = Ft. (14) 



An output signal from the oscillator such as X = cos0(t) will have a power spectrum 
consisting of a Lorentzian peak centered at Q (and, for more general signals, the harmonics) 

with 5*0 the spectral weight of the delta-function peak in the spectrum of the no-noise 
oscillator [19j. The width of the spectral peak is therefore equal to the phase noise strength 
F. Thus the tails of the spectrum away from the peaks decay as w"^; this is the white-noise 
component of the Leeson noise spectrum for oscillators [22] . Other noise spectra will lead to 
different power law tails. 

The relationship of the effective noise strength F acting on the collective phase G to 
the the strength of the noise / acting on each individual oscillator is given by projecting 
the individual noise components ^i(t) along the phase variable 9. Denoting the phases 9i 
by the vector 0, the tangent vector to the limit cycle (the zero-eigenvalue eigenvector) is 
given by eo = (1, 1, 1 ... , 1), choosing a convenient normalization so that the phase shifts 
corresponding to a time translation At are 66i = eo,jAt. Using the general results of refs.flSl- 
[20l or the simpler analysis for the present case sketched in Appendix the relationship 
is 

F = /, (16) 

with Bq^ the zero-eigenvalue adjoint eigenvector. Thus finding the broadening of the line due 
to the noise is reduced to calculating the adjoint eigenvector Bq . 

The Jacobean matrix J yielding the linear stability analysis of the phase dynamics about 
the fixed point phases O^'^^ defining the limit cycle is 

Jij = r'{6j'^^ — ol'^^) for i,j nearest neighbors, (17a) 

Ju = -J2r'i(^f-0h, (17b) 

nn j 

with other elements zero. The vectors eg, are defined by 

J • eo = 0, (17c) 
J"^ • eo"^ = 0, (17d) 

with j\j = Jji. 
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Note that I am treating noise perturbatively in the small noise limit and for a finite system: 
in this case the result is given by just the effect on the overall phase of the synchronized 
state, which is the zero-mode of the system. I am not considering modifications to the 
synchronized state due to the noise, such as changes in values of the critical disorder for 
synchronization, or changes in the nature of the synchronized state. In a finite system there 
will be barriers to such fiuctuations, and their rates will vary with the noise strength / as 
e~^/f with A some number depending on the states considered. These fiuctuations can 
therefore be ignored for small enough /. As the number of oscillators tends to infinity, some 
barriers will become very small, and the synchronized state may be significantly changed or 
even eliminated by the addition of noise [H], as for phase transitions in equilibrium systems 
at finite temperature. 



V. DISSIPATIVE COUPLING 

For purely dissipative coupling 7 = 0, the Jacobean J is symmetric and so the adjoint 
eigenvector is equal to the forward eigenvector which is the tangent vector defined by an 
infinitesimal time translation 



e, 



:o^ = eo = (l,l,l...,l). (18) 



This result is true for any antisymmetric coupling, and is not restricted to the nearest 
neighbor model. This immediately gives the result for the effective noise strength 

F = N-'f, (19) 

so that the frequency precision of the entrained state is enhanced by the factor N. Note 
that the enhancement does not depend on the degree of phase alignment quantified by the 

(s) 

magnitude of the order parameter \1/ = A^~^ ^ • e* ^ , which may be less than unity (i.e.. 



phases not fully aligned) even in the entrained state. The result Eq. (19) has been obtained 
previously 121 [3], although I believe the present derivation is more systematic, since it does 
not assume that the effect of the noise on the collective phase remains small at long times. 
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VI. DISSIPATIVE PLUS REACTIVE COUPLING 



For general coupling the Jacobean is not symmetric, and there is no obvious relationship 
between eg^ and eg in general. Physically, in situations where the entrained state consists 
of waves emanating from a source region of higher frequency oscillators, we might expect 
the frequency precision to be determined by fluctuations of only those oscillators in the core 
region that fix the frequency of the waves. This means that the reduction of the effective 
noise by averaging is only over this core region of oscillators, giving a poorer improvement 
of the frequency precision. I first demonstrate this result for a simpler "one-way" coupling 
function introduced by Blasius and Tonjes[T5] for which analytic solution is possible. I then 
derive the result for the general coupling function Eq. ^ assuming the disorder is small 
enough so that the phase difference between all nearest neighbor oscillators is small, in which 
case the approximation Eq. ^ may be used. The result depends on the mapping [1^ [15] 
of the solution for the entrained state onto the Anderson localization problem[23], and the 
known properties of the localized states in this problem|24t 1^. together with a relationship 
between and the localized states that I demonstrate. I also investigate one and two 
dimensional lattices numerically. 



A. One-way coupling 

Blasius and Tonjes[in] proposed a simple, exactly soluble model of a one dimensional 
lattice, with a nearest neighbor coupling function such that the phase of oscillator i is only 
influenced by its neighbors if their phases are ahead (all phase differences are assumed to be 
small so that this notion makes sense). I use the example 

f 7(e^<^ - 1) for > 0, 
m={ ) ^ ^ (20) 

I for < 0. 

The entrained solution is given hj Q = Um with m the index of the largest frequency in the 
lattice, and then the fixed point solution 6^'"^^ is constructed iteratively from oln' using 

q{s) ^ I C\ - 7"^ ln[l + l{(^m - UJ^)] for i > m, ^^^^ 
[ 6',^+! - 7""^ ln[l + 7(wm - uJi)] for i < m. 

The value chosen for 6m sets the overall phase B. An example of the entrained state for 
200 oscillators in a Id lattice is given in Fig. [T| showing waves emanating from the oscillator 
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FIG. 1. Entrained state for the one sided model Eq. (20) for 200 oscillators. Left panel: oscillator 



phase as a function of lattice site with 9m set to 0; right panel: gray scale plot of cos 9i as a function 
of time. The oscillator frequencies were taken from a uniform frequency distribution with width 
a = 0.5 and mean zero, and 7 = 1. 



with maximum frequency at m = 121. 

It is easy to see for this coupling function that the zero-eigenvalue adjoint eigenvector is 

^0,i ~ ^im-, (22) 

corresponding to the fact the phase 6^, is not coupled to either neighbor, since Om > Om±i- 
Thus the effective noise is given by F = /, and there is no improvement of the frequency 
precision, even though all the oscillators are entrained. This is because the single oscillator 
with maximum frequency determines the entrained frequency, and therefore the entrained 
frequency is as sensitive to noise as this single oscillator. 



B. General coupling 

I now consider the case of general coupling Eq. ^ in the limit of small enough disorder 
so that the small phase difference approximation Eq. (|6| approximation may be used. In 
this case, Blasius and Tonjes[T5] showed that the Cole-Hopf transformation 9i = 7~^lngj 
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maps the problem for the entrained state onto the hnear problem 

qi = Eqi = jUiQi + '^{qj - qi), (23) 

nn j 

with the eigenvalue E = 'yQ. This is equivalent to the tight binding model for a quantum 
particle on a random lattice, and the properties of the solution can be extracted in analogy 
with Anderson localization [23j. At long times the solution q(t) = qmaxgEnmxt corresponding 
to the largest eigenvalue E^ax will dominate. This gives the entrained state 

^f)=7-Mngr-, (24) 

with frequency fl = i?max/7- Anderson localization theory shows that q™'^^ may be chosen 
positive, and it has the form of an exponentially localized state centered on a region of the 
lattice with a concentration of larger frequency oscillators. The exponential localization of 
qmax corresponds to a roughy linear phase profile, again leading to waves propagating from 
a source, as shown in Fig. |2} 

I now analyze the frequency precision based on the properties of the solution q™^^ known 
from studies of the Anderson problem. 



1. One- dimensional lattice 

For a one- dimensional lattice with nearest neighbor coupling Eq. (|6|), the Jacobean matrix 



Eq.(|17|) for the stability analysis of the fixed point solution is the tridiagonal matrix with 



elements 

Ju±i = ef^^'^^-'^^'^ = gST/C'^ (25) 

Jii = ~Ji,i+l ~ Ji,i-1^ (26) 

except for the first and last rows corresponding to the end oscillators which only have one 
neighbor 

J 7 ^max /^max j j ^max /^max /n'7\ 

^^12 — —Jll — (I2 1^1 ^ 'JNN-1 — —'JNN — Qn-i/ In ' \-^') 

and all other elements zero. It is easily checked that (1,1,1...,1) is indeed the zero- 
eigenvalue eigenvector. The adjoint matrix has off diagonal elements 

A±i = QrVQT±i (28) 
11 
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FIG. 2. Entrained state for the general model Eq. ([5]) for a chain of 200 oscillators with nearest 
neighbor coupling, using the small phase difference approximation Eq. ([6]). Left panel: oscillator 
phase as a function of lattice site; right panel: gray scale plot of cos 9i as a function of time. The 
oscillator frequencies were taken from a uniform frequency distribution with width o" = 0.5 and 
mean zero, and 7 = 1. 



except for the first and last rows for which 

T t _ max / max t f _ max /„max /qqN 

and diagonal elements 

A^ = Jn. (30) 

with all other elements zero. The key result is that the (unnormalized) adjoint eigenvector 
can be found explicitly 

4 = (31) 

as can be confirmed by direct substitution. This simple result follows from the quotient 
form of the Jacobean matrix elements for the special form of the interaction Eq. (|6]). An 
example of the vector q^^^ and the adjoint eigenvector for the system of Fig. [2] is shown 
in Fig. [3j 



The noise reduction factor -F//, Eq. (16), is given by 



12 



21 2- 



(32) 




80 



100 120 
Index 



140 



o 
^— ' 

o 
c 

'-t—l 

o 

13 

■D 


i_ 


CO 




0.1 r 



YO 



FIG. 3. Left panel: Localized solution q™^^ (blue circles) and zero-eigenvalue adjoint eigenvector 
(black squares) for the system of Fig. ^ (The normalizations are chosen for the plot so that 
the largest element of each vector is 1.) Only the portion of the L = 200 lattice where the elements 
have appreciable size is shown. Right panel: scaling of the noise reduction factor F/f deduced 



from Eq. (32) with 7cr. Each point is the average of 1000 realizations of the random lattice of 
frequencies. 



This equation directly relates the improvement in frequency precision to the solution of 



the linear Anderson problem Eq. (23). The expression Eq. (32) for the noise reduction is 



the inverse participation ratio p^^ of the vector q™^^ of the linear localization problem. 
This can be used to define the radius of the localized state r = p/2. Thus I find that the 
noise reduction factor is given by the size of the source of the waves, rather than by the 
total number of oscillators, giving an improvement in frequency stability that is significantly 
worse for a large number of entrained oscillators. The size of the source is defined precisely in 
terms of the participation ratio of the maximum energy localized state of the corresponding 
Anderson problem. For the system in Figs. 2][3, p ~ 9.36. In this example, the frequency 
precision would not be improved by increasing the number of oscillators beyond about ten. 



From Eq. (23) it is clear that the noise reduction factor F/f depends on the parameters 



of the model only through the product 70", for the approximation to r(0) used. Within this 



approximation, the scaling found from numerical solutions of Eq. (23) for a one dimensional 
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lattice is shown in Fig. |3| The calculations were done for chains of length 100-1000, and 
the results were insensitive to the length providing it is much larger than the width of the 
localized state. A power law F// oc (7cr)°'^ is a good fit to the calculated results over the 
range considered. 

2. d- dimensional lattice 
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FIG. 4. Source and waves for a 60 x 60 two dimensional lattice of oscillators with 7 = l,cj = 2. 
Left panel: zero-eigenvalue adjoint eigenvector showing the effective size of the source of the 

(s) 

waves in the entrained state; right panel: grey scale plot of cos 91 J giving a snapshot of the waves 
emanating from the source. 



The same argument applies to general dimension, although the structure of the Jacobean 
matrix is no longer tridiagonal. Choose any convenient labeling of the oscillators 9i,i = 
1 . . . N. Nearest neighbor oscillators will not in general be adjacent in the list. However the 
Jacobean and its adjoint are still defined by 



Jij = q. 



max /„max 
3 



T t max / max 

•J'^ — Hi iHj 



= J- 



for ij nearest neighbors 
for ij nearest neighbors 



(33) 
(34) 
(35) 



nn J 



with other elements zero. The eigenvectors eo, are as before, and the expression Eq. (32) 
for F in terms of the inverse participation ratio is unchanged. Thus I expect the noise 



14 



reduction factor to scale as with r ~ p^^'^ the radius of the maximum energy locahzed 
state in the d-dimensional Anderson locahzation problem. 

Figure |4] shows an example of a target wave entrained state for a 60 x 60 two dimensional 
lattice. The left panel is the adjoint eigenvector defining the source: the inverse participation 
ratio, yielding the improvement in the frequency precision, is 54.1 (cf. = 3600). The right 
panel is a plot of cos6lj\ calculated from the Cole-Hopf transformation and the eigenvector 
qmax^ giving a snapshot of the waves in the entrained state. 



C. More general systems 



The results Eq. (16) and Eq. (31) remain valid for a more general coupling yielding the 



equations for the phase dynamics 

ei = coi + J2 - e,) + ut), i = i,...,n. (sg) 

with Kij a symmetric matrix giving the strength of the coupling between oscillators i and 
j. The small phase difference condition so that Eq. ^ may be used is now that the phase 
difference between any two oscillators with nonzero K^j be small in the entrained state. The 



same analysis leads to the relationship Eq. (32) between the noise reduction factor and the 



participation ratio of the largest E eigenvector q™**^ of the corresponding linear problem 

Eqi = jUiqi + ^ Kij{qj - qi). (37) 

Note that although the strength of the coupling can be different for different paris of oscil- 
lators, the forvfi of the coupling Eq. (|5|, and in particular the ratio of reactive to dissipative 
components, must be the same for this simple analysis to apply. 



One generalization Eq. (36) allows is to a lattice of oscillators with short range, but 
not just nearest neighbor, interactions. The scaling of the noise reduction with 70" will be 
the same as for nearest neighbor interactions, since the scaling properties of the Anderson 
problem are the same for these two cases. More generally, the method reduces the problem 
of calculating the improved frequency precision in the entrained state of a complex network 



of oscillators to solving the linear problem Eq. (37) for the network architecture and coupling 
parameters K^. 
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VII. DISCUSSION 



The major result of this paper is that for oscillators on a lattice with short range coupling 
including a reactive component, the improvement of the frequency precision due to synchro- 
nization is limited to a factor given by the number of oscillators in the core source region of 
the waves that form the entrained state, rather than a factor equal to the total number of 
oscillators, as is the case for purely dissipative coupling. I showed this result explicitly for 
the phase reduction description, in the limit of small enough disorder, or strong enough cou- 
pling, so that the phase differences between interacting oscillators are small in the entrained 
state. The size of the core region is given by the extent of the localized ground state of the 
corresponding linear Anderson problem, onto which the nonlinear phase equation is mapped 



by a Cole-Hopf transformation. The precise relationship is Eq. (32) relating the reduction 
in the phase noise to the inverse participation ratio of the localized state. This relationship 
remains true for general networks of oscillators providing the small phase difference approx- 
imation Eq. ^ applies for all interacting pairs of oscillators, and reduces the calculation of 
the frequency precision to the corresponding linear Anderson problem on the network. 



Within the small phase difference approximation, the entrained state of waves propagating 
from the localized source is the unique state at long times. However, for the phase equations 
with the full coupling function Eq. ([s]) other states may result depending on the initial 
conditions. This is particularly evident for two dimensional lattices, where spiral states 
are seen in numerical simulations starting from particular initial conditions |14j. Due to 
the topological constraint of integral 27r phase winding around the center, such a structure 
survives at long times, unless the core migrates to an open boundary. A second consequence 
of the topological structure is that there are necessarily large phase differences between 
nearest neighbor phases in the core, so that the small phase difference approximation breaks 
down. In the spiral state all the oscillators are again entrained to a single frequency, that 
probably depends just on the oscillators in the core region. It would be interesting to extend 
the analysis of the present paper to these spiral states. 
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Appendix A: Derivation of the phase equation 



In this appendix I give a brief derivation of the stochastic phase equation Eq. (11). This 
equation follows from the general results for limit cycles of refs.[T2]and|2Dl but the derivation 
is simpler for the phase reduction description. I start from Eq. ([T]) 

e, = uj, + Y. - + vUt), t = l,...,N, (Al) 

nn j 

introducing a perturbation parameter r] to label the small noise. I will develop a perturbation 
expansion in rj, and set r/ — )■ 1 at the end of the calculation. I will expand to first order 
in 7] to extract the phase diffusion: continuing to second order would be needed to find the 
Lamb shift like correction to the mean frequency. 

At zeroth order in rj the solution is the no-noise solution Eqs. (|8 10 ), with Go an arbitrary 
constant. In the presence of order t] noise I expect this overall phase to evolve on a slow 
time scale T = rjt, Oq = 6o(T), so that the phases will be given up to order r] by 

9i (t) = el'^ + ^]t + Bo (T) + (t) + . . . , ( A2) 

with 6^^\t) a correction to be found. Expanding the equation of motion up to order t] the 
equation for this correction is 

^f^-E^^^-^f = -(©oeo.-e.(t)), (A3) 



with Jij the Jacobean Eq. (17), Q'q = dQo/dT, and eo = (l,l,l...l). Components of 
6^^^ along eigenvectors of J with negative eigenvalues will have some finite value given by 
inverting this equation. However, for the component along the zero-eigenvalue eigenvector, 
there is no restoring force, and any nonzero value of the right hand side will lead to large 
values of 6^^^ at large times, violating the assumption that 6'^^^ gives a small correction. 
This component is extracted by multiplying on the left by the adjoint eigenvector e^^ since 
Bq ■ 3 = 0. This leads to the solvability condition for 6^^'' to remain finite 



t 

=0 ■ ^0 



0'o = 



Returning to the original variables and setting r/ — )■ 1 gives the stochastic equation for the 
overall phase 

Q{t) = n + E{t), (A5) 
17 



with the correlation function of the effective noise 



E.,4eo|,fe(t)e.(t')) 



■ eoj 



(A6) 



giving Eq. (12) with Eq. (16) for equal, uncorrelated white noise of strength / acting on each 



individual oscillator. Note that the result does not depend on the choice of normalization 



for Bq^. a specific normalization choice for eg was made in setting up Eq. (A3). 
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